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ABSTRACT 



Strange baryons have long been known to exhibit a leading particle effect. 
A recent comparison of S~ production in tt~ , n, and E~ interactions with 
nuclei show this effect clearly. These data are supplemented by earlier mea- 
surements of and Q production by a beam. We calculate the S~ 
and fl xp distributions and nuclear dependence in hA interactions using the 
intrinsic model. 
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1 Introduction 



Leading particle effects, flavor correlations between the flnal-state hadron and the pro- 
jectile valence quarks, have long been observed in strange particle production. Although 

many experiments have recently focused on leading charm production [1. 2, 3, 4, 5, 
6, 7, 8, 9, 10], the first data involved strange particles [11, 12, 13, 14, 15, 16]. With 
new data from the WA89 collaboration on S~ production by 7t~ , n, and S~ projectiles 
on nuclear targets [17], in addition to E~ production data from S~ beams [14], doubly 
strange hadron production can be studied as a function of the number of strange valence 
quarks in the projectile. We compare our model calculations to both the xp distribu- 
tions and the integrated A dependence reported by WA89 [17]. We also discuss H~ and 
fl production by the S~ beam [14]. 

The WA89 collaboration used carbon, C, and copper, Cu, targets to study the A 
dependence of 'E.~{dss) production by 7r~{ud), n{udd), and H~{dds) beams [17]. The 
negative beams, vr^ and E^, had an average momentum of 345 GeV with a 9% mo- 
mentum spread. The neutron beam had a lower momentum with a larger spread than 
the negative beams — the average momentum was 260 GeV with a 15% variation. The 
detected E~ was in the forward xp region, xp > 0.05, with low transverse momentum, 
Pt < 2.5 GeV. The data were parameterized in the form 

"^"^ oc (1 - xfTc-^'^ . (1) 



dpji^dxp 

The pion and neutron results agree with the functional form of Eq. (1) over all xp. For 
the pion, a = 3.8 ± 0.3 for C and 4.1 ± 0.3 for Cu while for the neutron a = 5.0 ± 0.3 for 
C and 4.8 ± 0.3 for Cu. These results are consistent with expectations from spectator 
counting rules [18], da/dxp oc (1 — a;^)^"""^. With an incident gluon, = 2 for pions 
and 3 for neutrons, consistent with no leading particle effect for projectiles with zero 
strangeness. There is no significant A dependence of the exponent a. 

On the other hand, the data cannot be fit to Eq. (1) for xp < 0.4. In the large 
Xp region, a = 2.08 ± 0.04 for C and 1.97 ± 0.04 for Cu. These results indicate a very 
hard xp distribution, inconsistent with the counting rules even for a valence quark since 
Us — 2 gives {1 — xp)^. In addition, at xp < 0.4, the distribution is independent of 
Xp for both targets. Thus these data show a strong leading particle effect since the S~ 
has two valence quarks in common with the S~. The statistics are also sufficient for an 
observable A dependence in the fitted values of a. 

The integrated A dependence was also reported by WA89 [17]. The A dependence 
of the total cross section is often parameterized as 

(7pA = (^ppA"' . (2) 

The integrated a for E~ production of S", a = 0.679 ± 0.011 [17], is in relatively good 
agreement with previous fits. However, the pion and neutron data show a closer-to-linear 
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A dependence, a = 0.891 ± 0.034 and 0.931 ± 0.046 respectively. WA89 attributes this 

difference to the fact that two ss pairs must be produced to make the final-state S~ and 
ss pair production would be suppressed relative to light qq production. 

WA89 has also measured the dependence oi a on xp. This dependence, a{xF), was 
previously reported for a wide range of hadron projectiles [19]. For non-strange hadrons 

and hadrons with a single strange quark, there is a common trend with xp. At Xp = 0, 
a ~ 0.8 and decreases to ^ 0.5 at large overall decrease of ~ A^^^ for < xp < 1. 

The S*^, the only doubly-strange hadron included in Ref. [19], is an exception. In pA 
interactions, the S° has a larger value of a at low xp [16]. A similar effect is observed for 
S~ production by WA89. Their measurements of a{xp) for S~ from pion and neutron 
beams show that a 1 for xp 0.05, decreasing to a ~ 0.7 at higher xp. Thus the 
decrease of a with xp is also A^^^ in this case although the actual values of a arc larger 
than those for lighter hadrons [19]. However, for S~-induced S" production, a ~ 0.7 
almost independent oi xp. 

The other data we consider are S~ Be — > S~, fi(sss) at 116 GeV, measured by Biagi 
et al. [14]. In this case, the final-state 'E~ xp distribution increases with xp^ as does 
the VL Xp distribution. This increase could be due in part to the use of an invariant 
parameterization [14], 

oc (1 - xpT'e-^'^r , (3) 

which fits the S~ data at xp > 0.5 but only approximately fits the fl data in this limited 
region. The exponent a' was fit in two pf. intervals, p'^ < 0.4 GeV^ and 0.4 < p"^ < 2.9 
GeV^, yielding a' = -0.45 ± 0.02 and -0.18 ± 0.03 respectively. Between the most 
central measurement, xp = 0.15, and the projectile fragmentation region, xp — 0.85, 
the E~ cross section increases by a factor of ~ 40 in the low p'^ interval. 

A comparison of these results with incident proton data [14], pA — > E~X [11, 20] and 
pA — > pX [21], showed that, at low xp, the S~ production cross section is essentially 
independent of the projectile while, at high xp, the S~ and p scattering cross sections 
are similar. This behavior supports valence quark domination at high xp. The structure 
of the Q Xp distribution is similar: it is of the same order of magnitude as pA — > flX 
[11] at low Xp but is similar to singly strange baryon production by protons, pA — > AX 
and E+X, [11, 20, 22] at high xp. 

Since only one target was used, a = 0.6 was assumed in Eq. (2) to obtain the per 
nucleon cross sections. This extrapolated cross section is a factor of 1.5 — 2 higher than 
those on hydrogen targets [14]. An extrapolation with a — 1 gives better agreement 
with the hydrogen target data, at least for E~ production. 

We employ the intrinsic model [23, 24, 25, 26, 27] as developed for strangeness pro- 

^For the two parameterizations to be equivalent, the right-hand side of Eq. (1) should be multiplied 
by 2El^ to obtain the invariant cross section. 
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duction in Ref. [28]. In the intrinsic model, a hadron can fluctuate into Fock state 
configurations with a combination of light and strange quark pairs. The heavier quarks 
in the configuration are comoving with the other partons in the Fock state and thus 
can coalesce with these comoving partons to produce strange hadrons at large xp- The 
model combines leading-twist production of ss pairs with intrinsic Fock states with up 
to nine particles. Thus coalescence production of the Q from a proton is possible. 



2 Leading-Twist Production 



We calculate leading- twist strangeness in perturbative QCD, assuming the strange quark 
is massive. When the projectile has nonzero strangeness, we also consider the possibility 
of flavor excitation. We choose proton parton distribution functions with the lowest pos- 
sible initial scale ^Uq so that ml > /Xq. Therefore the baryon parton distribution functions 
are based on the GRV 94 LO proton parton distributions [29]. We use the most recent 
pion parton densities by Gliick, Reya and Schienbein [30]. To be conservative, we as- 
sume that the scale /i at which the strong coupling constant and the parton densities are 
evaluated is fj, — 2mT where mr = \lv\ + "^s — 500 MeV. The xp distributions, 

obtained by integrating Eq. (1) or (3) over p^, are dominated by low pt production. 

We treat the strange quark as heavy, as in Ref. [28], rather than as a massless 
parton in jet-like processes. Treating the strange quark as a jet means selecting a 
minimum to keep the cross section finite. A large minimum p^ compatible with hard 
scattering is incompatible with the assumption of intrinsic production, inherently a low 
Pt process [25]. Strange hadrons can either be produced directly in jet production or 
by the fragmentation of light quark and gluon jets. However, there is no indication that 
these data originate from jets. 

The leading-twist xp distribution of heavy quark production [31] is denoted by F, 

F = = — / dZ3dy2dpTXaXbHAB{Xa:Xb,ljL )— '- , (4) 

aXp Z J Hii Z3 

where A and B are the initial hadrons, a and b are the interacting partons, 1 and 2 are 
the produced strange quarks and 3 is the final-state strange hadron S. The xp oi the 
detected quark is xp — 2mr sinh y/y/s where y is the rapidity of the quark and ^/s is the 
hadron-hadron center of mass energy. We assume the simplest possible fragmentation 
function, 

Ds/s(z)^Bs5{l-z) , (5) 

with Bs = 0.1, assuming that all 10 ground-state strange hadrons are produced at the 
same rate to leading twist [28]. This choice of Ds/s gives the hardest possible leading 
twist Xp distribution. 
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The convolution of the subprocess qq annihilation and gluon fusion cross sections 
with the parton densities is included in HAsixa, ^b, l-i'^), 

where q — u, d, and s. Although including the s quark in the sum over q in Eq. (6) 
could lead to some over counting, the strange quark contribution to F from non-strange 
projectiles is negligible, less than 0.01% for neutron and pion projectiles. It is some- 
what larger for strange projectiles, 2.5% for the S~ and 5.6% for the S~ but it is only- 
significant at large xp- 

Hyperon parton distributions can be inferred from the proton distributions [27] by 
simple counting rules. We can relate the valence s distribution of the S^, /Jj" , to the 
proton valence d distribution, /^^, and the valence d distribution in the E~, , to the 
valence u in the proton, /^^, so that 

f^dxfZ{x,i/) = J^'dxfl{x,^i')^l, (7) 

£dxfl-{x,i,') = £dxfl{x,i,')^2. (8) 

We also identify the up quark in the sea of the E~ with the strange quark in the 
proton sea, /J {x, /j^) — ff{x, //^). Similar relations hold for the antiquark distributions. 
Likewise, for the S~, we relate the valence s, , to the valence u in the proton, , 
and equate the valence d distributions in both baryons so that, 

£dxff;{x,^') = £dxfl{x,^') = 2, (9) 

£dxfl-{x,i,') = jyxflix,!,')^!. (10) 

Here also, (a^, /^^) = fsi^ii^'^)- The gluon distributions are assumed to be identical 
for all baryons, fg=ff ^ fg ■ '^^^ leading order subprocess cross sections for 
heavy quark production can be found in Ref. [32]. The fractional momenta carried by 
the projectile and target partons, Xa and Xb, are Xa — {ttit/ y/s){e^ + e^^) and x^ — 
{mT/^/^){e-y + e-y^). 

We have assumed only gg and qq contributions to massive strange quark production. 
We have also checked how the xp distribution would change if the strange quark was 
treated as massless and all 2 — > 2 scattering channels were included. Jet production of 
s quarks is through processes such as g^s — > gs, qs qs and qs qs. (Similarly for the 
s.) Including these processes increases the cross section by a factor of 4 — 8. While this 
factor is not constant, it increases rather slowly with a;i? so that the difference in shape 
is only important in the region where intrinsic production dominates, as discussed later. 
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Contributions from massless 2 — > 2 scattering increase more rapidly at Xi? > for 
strange projectiles because the contribution from, for example, {xa)fg{xb)i domi- 
nates the scattering cross section. In the infinite momentum frame, ff^ = f^^, see 
Eq. (8), and ff^ is large at large Xa while increases as Xb decreases. To take this 
into account quantitatively, we have incorporated "flavor excitation" of massive strange 
valence quarks. The excitation matrix elements for massive quarks are found in Ref. [33]. 
The flavor excitation cross section has a pole when — > so that a cutoff, PT^^^^ is re- 
quired to keep this cross section finite, as in jet production. We employ PT^in = = 1 
GeV. The leading-twist fusion cross section for strange projectiles is then augmented by 

where 

H'ABKA.t^')-ft{<^t^'){Uf!{^t^^ > (12) 

x'^ = {rriTey + Pt^^'^) / y/s and x[ = {mxe~y + pxe^"^'^) / y/s ■ Note that there is no overlap 
between the processes included in Eqs. (12) and (6) and thus no double counting. This 
excitation mechanism is effective only for hadrons with a strange quark in the final state 
and thus does not affect the distributions with a produced s. 

To summarize, for strange and antistrange final states produced by non-strange 
hadrons, 

-j^^-^ = F , 13 
dxp dxp 

as in Eq. (4). This relation also holds for antistrange final states from strange hadrons. 
However, for strange hadron production by hadrons with nonzero strangeness, we also 
consider 

^^^j^^^^p j^X (14) 
dxp dxp dxp 

where flavor excitation, Eq. (11), may play a role. 

We remark that the role of flavor excitation in heavy quark production, as outlined 
in Ref. [33], is questionable. It was flrst proposed as a leading order contribution to the 
total cross section and, as such, could be rather large if the heavy quark distribution in 
the proton is signiflcant. However, the proton heavy quark distribution is only nonzero 
above the threshold mq. In addition, parton distribution functions are deflned in the 
inflnite momentum frame where the partons are treated as massless. Later studies at 
next-to-leading order (NLO) [34] showed that these excitation diagrams are a subset 
of the NLO cross section and suppressed relative to fusion production. They are only 
a small fraction of the heavy flavor cross section and thus play no signiflcant role at 
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Figure 1: Leading-twist strange quark production in (a) n^p interactions at 345 GeV, (b) 
np interactions at 260 GeV, (c) T,~p interactions at 345 GeV, and (d) E~p interactions 
at 116 GeV. The solid curves are the fusion results, F. For projectiles with valence 
strange quarks, the excitation contributions, Xp^ with pr^in ~ 1 GeV, are shown in 
the dashed curves. The dot-dashed curves are the total, F + Xp . 

min 
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low energies. Strange hadron production at large xp is then an important test of the 
excitation process. 



The leading- twist xp distributions with all four projectiles are shown in Fig. 1. We 
compare the fusion xp distributions, F, in 7r~p and interactions at 345 GeV, np 
interactions at 260 GeV, and S~p interactions at 116 GeV, corresponding to the energies 
wc investigate. The 7r~p distribution is the broadest because the f^^ {^a)fu^{xb) contri- 
bution hardens the xp distribution at large xp where the qq channel dominates. When 
excitation is considered, as in H^p and E^p interactions, the xp distribution is hard- 
ened, particularly through the {^'a)fg{^b) /f {x'^)fg{x[) contributions. These 
dominate at large xp where fs{x'g) is large for valence strange quarks and f^{x'^) is 
large at small x[. The effect is even stronger for the S~ since it has two valence strange 
quarks. Note that, except at small xp, the total leading-twist cross section is equivalent 

to Xr,„ 

" min 

We can obtain an approximate estimate of the exponent a from Eq. (1) from the 
average xp, (xp), where 



When averaged over xp > 0, the values of a obtained are larger than those measured 
by WA89, as expected. For the pion and neutron beams, a — 5.2 and 8.7 respectively. 
Strangeness production by strange hadrons including fusion alone also gives large values 
of a, 9 for the E~ and 7.4 for the S~. The xp distribution of strange quarks produced 
by flavor excitation is considerably harder, a = 3.2 for the S~ and 2.2 for the S~. 
Combining the two contributions, as in the dot-dashed curves in Fig. 1(c) and (d), gives 
a somewhat larger value of a than for flavor excitation alone, a — 4.7 and 2.9 for E~ 
and S~ beams respectively. The values of a obtained from Eq. (15) are all much higher 
than those obtained from the data. Thus the leading twist results alone cannot explain 
the shape of the measured E~ xp distributions. 



3 The Intrinsic Model for Strangeness 

We now briefly discuss the intrinsic model for strangeness production, described in detail 
for n~p interactions in Rcf. [28]. Since all the data is at > 0, we only discuss intrinsic 
production in the projectile. 

The hadron wavefunction is a superposition of Fock state fluctuations in which the 
hadron contains one or more "intrinsic" QQ pairs. These pairs can hadronize when the 
hadron interacts, breaking the coherence of the state. The model, first developed for 
charm [23, 24], gives heavy quarks a larger fraction of the projectile momentum due 
to their greater mass. The strange quark is lighter so that the momentum gained is 



1 



(15) 
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not as large. However, the intrinsic strangeness probability is larger, ~ 2%. For 
simplicity, we assume that the intrinsic probabihties are independent of the valence 
quark content of the projectile. Then is identical for nucleons and hyperons. The 
Fock state probabilities for up to 3QQ pairs where at least one QQ pair is strange are 
given in Ref. [28]. 

In this paper, we focus on E~{dss) and fl{sss) production from n~(ud), n{udd), 
J]~{dds) and projectiles. The produced S~ shares one or more valence quarks 

with the projectile. We study fl production only by S~ projectiles, with two valence 
quarks in common. 

Once the coherence of the Fock state is broken, the partons in the state can hadronize 
in two ways. The first, uncorrelated fragmentation of the strange quark, is the same basic 
mechanism as at leading twist. However, when the Fock state fluctuation includes all 
the valence quarks of the final-state hadron, these quarks, in close spatial proximity, can 
coalesce into the final-state hadron and come on shell. Thus, to calculate the full strange 
and antistrange hadron xp distributions in the intrinsic model, we include uncorrelated 
fragmentation of the strange quark in every state considered and coalescence from those 
states where it is possible. In Ref. [28], a comparison of the model with strange baryon 
asymmetries suggested that fragmentation may not be an effective mechanism because 
when the Fock state has minimal invariant mass, fragmentation may cost too much 
energy. This conclusion needs to be checked against inclusive xp distributions over a 
broader xp range. 

In principle, the parton distributions of the hadron can be defined through such a 
Fock-state expansion [35]. In each fiuctuation, only the mass distinguishes the fight and 
heavy quark distributions. Thus it is not really possible in a given state to separate the 
"valence" and "sea" distributions. All are similar as long as the quarks are light. One 
distinguishing feature is our assumption that only strange quarks can produce strange 
final-state hadrons by uncorrelated fragmentation. Thus with hyperon projectiles, un- 
correlated fragmentation may also be possible from Fock states with only fight QQ pairs. 
These intrinsic light quark states must be included in the total probability, as described 
in Ref. [28] . The probabilities for these states must also be defined. We assume 



We further assume that the probabilities for the meson Fock configurations are equal to 
the baryon probabilities. 

We have only taken the 10 strange ground state hadrons and antihadrons into ac- 




(16) 



(17) 



(18) 
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count. We assume that each hadron has a 10% production probabihty from fragmenta- 
tion, neglecting the particle masses. The final-state xp distribution is then equivalent 
to that of the s or s quark. For coalescence, we count the number of possible strange 
and antistrange hadron combinations that can be obtained from a given Fock state and 
assign each strange hadron or antihadron that fraction of the total. The possible num- 
ber of strange hadrons is greater than the number of possible strange antihadrons. We 
clearly err in the overall normalization by simply including the ground state strange 
particles. However, the higher-lying resonances have the same quark content with the 
same fragmentation and coalescence distributions since all properties of the final-state 
hadrons except their quark content are neglected. 

To obtain the total probability of each strange hadron to be produced from projectile 
hadron, h, in the intrinsic model, we sum the probabilities over all the states with up to 
3 intrinsic QQ pairs. Thus 

^^S _y^y^y^y^/o/ ^ i(mss)(m„u)(mdd) i(mss)(m„u)(mdd) I 

dxF rit^tr VIO dxF dxF I - ^ ' 

To conserve probability, (3 — 1 when the hadron is only produced by uncorrelated 

fragmentation and 0.5 when both fragmentation and coalescence are possible. When we 
assume coalescence production only, we set P"^ = and j3 = 1. The weight of each 
state produced by coalescence is ^ where ^ = Q when S is not produced by coalescence 
in state |n^ms(ss)m„(MM)md(dd)). The number of up, down and strange QQ pairs is 
indicated by m„, and respectively. The total, -|- -|- = m, is defined as 
m = {n — n^)/2 because each Q in an n-parton state is accompanied by a Q. For baryon 
projectiles, n = 5, 7, and 9 while for mesons n = 4, 6, and 8. Depending on the value of 
n, rui can be 0, 1, 2 or 3, e.g. in a \uudssdddd) state, = 0, = 2 and = 1 with 
m — 2). Note that = is only possible when h is strange since no additional ss pairs 
are thus needed to produce some strange hadrons by coalescence. The total probability 
distributions, dPg/dxp, for S — E~ and ft are given in the Appendix. 



4 A Dependence of Combined Model 



The total xp distribution for final-state strange hadron S is the sum of the leading-twist 
fusion and intrinsic strangeness components, 

dxp dxp dxp 

The leading-twist distributions are defined in Eqs. (13) and (14). The total intrinsic 
cross section, dafq/dxp, is related to dPg/dxp by 

^^iQ _ in f^' dPi 

dxp mi dxp • ^^^^ 
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We assume that the relative A dependence for leading-twist and intrinsic production 
is the same as that for heavy quarks and quarkonia [25, 27, 31, 36]. The A dependence 
of the two component model is 

ahA = Af^ait + A^aiQ (22) 

where the combination of the two terms should approximate A" in Eq. (2). There are 
no strong nuclear effects on open charm at leading twist so that the A dependence is 
linear at a;^ ~ 0, a = 1, [37], dropping to a = 0.77 for pions and 0.71 for protons [31, 38] 
at higher xp where the intrinsic model begins to dominate. We assume that P — 1 and 
7 = 0.77 for pions and 0.71 for all baryons. Thus, 

A^-^ ^ ^-V3 as XF^l. (23) 

This relative A dependence, similar to that discussed earher for hght hadrons [19], is 
included in our calculations. The proton and neutron numbers are taken into account in 
the calculation of the leading-twist cross section. This isospin effect is small for fusion, 
F, which is dominated by the gg channel. In perturbative QCD, (3=1 could be modified 
by changes in the nuclear parton distributions relative to the proton [39]. However, the 
scale for our perturbative calculation is too low for such models of these modifications 
to apply [40, 41] and are not included in our calculations. 



5 Results 



We begin by comparing the model to the WA89 pion data in Fig. 2. These data do 
not strongly distinguish between leading-twist fusion and the full model. The intrinsic 
results do not significantly depend on uncorrelated fragmentation. All three curves agree 
rather well with the data, primarily because the fusion xp distribution is already fairly 
hard. Then the intrinsic contribution is a small effect even though the d valence quark is 
common between the 7r~ and the S~. This is perhaps due to the fact that P^^ is already 
rather small, ~ 0.6%. We note that the calculated total cross sections agree with the 
measured cross sections to within better than a factor of two despite the rather large 
uncertainties in the calculations. 

Even though the intrinsic contribution is relatively small, it significantly affects the 
value of a obtained from Eq. (15). The difference between the a values found without 
and with uncorrelated fragmentation in the intrinsic model is negligible for the pion 
beam. We find a ^ 4.1 for the C target and 4.3 for the Cu target relative to a = 5.2 for 
leading twist alone. These results are within the errors of the WA89 fit to their data. 
The agreement is especially good since the two-component model does not give a smooth 
falloff as a function oi xp that can be easily quantified by a single exponent. 

It is also possible to calculate a{xp) and the x^-integrated a from the ratios of the 
distributions. The calculations including both uncorrelated fragmentation and coales- 
cence generally give a smaller value of a and, hence, a stronger A dependence. This 
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Figure 2: The model is compared to the 345 GeV WA89 pion data on (a) C and (b) Cu 
targets. The dotted curves are leading-twist fusion, F, alone, the dashed curves include 
uncorrelated fragmentation and coalescence, and the solid curves include coalescence 
alone. The data sets have been normalized to the cross section per nucleon. The curves 
are normalized to the data at xp — 0.15. 



is because fragmentation peaks at low xp, influencing the A dependence sooner than 
coalescence alone which is only signiflcant at intermediate xp. Thus a{xp) ~ 0.9 for 
fragmentation and coalescence while a decreases from ~ 1 at low xp to 0.86 at high xp 
with coalescence alone. The integrated values are 0.93 and 0.98 respectively, somewhat 
higher than the WA89 result but with the same general trend. 

The overall agreement with the total cross section is not as good for the nA data, 
shown in Fig. 3. Surprisingly, the distribution including uncorrelated fragmentation 
agrees best with the data. This is perhaps because the energy of the secondary neutron 
beam is least well determined. The energy spread is 15% compared to 9% for the pion 
and S~ beams. A small energy variation can have a large effect on the leading-twist 
cross section. A 15% increase in the average neutron energy, from 260 GeV to 300 GeV, 
increases da^/dxp by 80% at Xi? ~ 0.25 while the intrinsic cross section is essentially 
unaffected. Such a shift in the relative leading-twist and intrinsic production rates could 
easily reduce the effect of coalescence alone to be more compatible with the data. The 
uncertainty in the energy of the pion beam has a much weaker effect on the flnal result 
because the intrinsic contribution is already small, as is obvious from Fig. 2. 

The calculated exponents a are larger for the neutron than the pion, in agreement 
with the WA89 measurements [17]. We find a ^ 4.9 for C and 5.8 for Cu. Typically the 
value of a obtained for fragmentation and coalescence is larger than that for coalescence 
alone since eliminating the fragmentation contribution tends to increase (xp). The 
stronger A dependence assumed for the intrinsic model has the effect of increasing a for 
larger nuclei. Thus the a found for the carbon target agrees rather well with the WA89 
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Figure 3: The model is compared to the 260 GeV WA89 neutron data on (a) C and (b) Cu 
targets. The dotted curves are leading-twist fusion, F, alone, the dashed curves include 
uncorrelated fragmentation and coalescence, and the solid curves include coalescence 
alone. The data sets have been normalized to the cross section per nucleon. The curves 
are normalized to the data at xp — 0.15. 

data while the copper data suggest a harder distribution than our calculation implies. 
There is, however, a stronger A dependence in the fallofT with xp than in the data 
which cannot distinguish between the values of a determined for the two targets. This 
stronger dependence is reflected in the values of a obtained, 0.87 when fragmentation 
and coalescence are included and 0.95 with coalescence alone. 

Even though there is some A dependence in the model calculations, the relatively 
small intrinsic contribution to the pion and neutron data leads to a rather weak overall 
A dependence. Dominance of the leading- twist cross section at low to intermediate xp 
results in a nearer-to-linear integrated a, as observed by WA89 [17]. 

We now turn to production by the E~ where the A dependence can be expected 
to be stronger. For the first time, we have a valence strange quark in the projectile so 
that we can compare the effectiveness of fusion alone with flavor excitation. We can also 
better test the importance of uncorrelated fragmentation because coalescence production 
is already possible in the 5-parton Fock state, \ddsss). 

Our results are collected in Fig. 4. We first discuss the importance of uncorrelated 
fragmentation to leading- twist fusion, F, alone. Fig. 4(a) and (b). The leading-twist 

contribution is rather steeply falling. Including both uncorrelated fragmentation and 
coalescence broadens the xp distribution but cannot match the hardness of the measured 
Xp distribution. Eliminating the fragmentation contribution produces a much harder 
distribution for xp > 0.15, matching the shape of the data relatively well. 
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Figure 4: The model is compared to the 345 GeV WA89 data on C and Cu targets. 
In (a) and (b), the leading twist contribution is F while in (c) and (d), flavor excitation 
is also included, F + X„„ . The dotted curves are for leading-twist alone, the dashed 

' ^ min 

curves include uncorrelated fragmentation and coalescence and the solid curves include 
coalescence alone. The data sets have been normahzed to the cross section per nucleon. 
The curves are normalized to the data at xp — 0.15. 
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The calculations are all normalized to the xp — 0.15 point to better compare the 
shapes of the distributions. Fragmentation gives better agreement at low Xp because 
uncorrelated fragmentation peaks close to xp ~ 0, filling in the low to intermediate 
xp range. Coalescence, on the other hand, always produces strange hadrons with 
{xp) > 0.3, broadening the distribution only in this region. Thus without fragmen- 
tation the calculations overestimate the data Bit xp ~ 0. The data seem to indicate 
that uncorrelated fragmentation is not an effective mechanism for intrinsic production, 
in agreement with the conclusions of Ref . [28] . 

The agreement with the solid curves in Fig. 4(a) and (b) is good but not perfect. 
The calculation overestimates the data at high xp. Recall that for the neutron, the 15% 
spread in the beam momentum could result in an overestimate of the intrinsic contribu- 
tion, as previously discussed. Although the possible spread in the S~ beam momentum 
is smaller, it could affect the relative intrinsic contribution at low to intermediate xp. At 
large xp, the effect on the shape would be negligible because the intrinsic contribution 
dominates. Thus, given the inherent uncertainties in the model and in the data, the 
agreement is rather satisfactory. 

We have obtained the value of the exponent a from (xp), both over all xp and for 
Xp > 0.1, avoiding the strong change in slope of the solid curves when coalescence alone 
is included in the intrinsic result. When the entire forward xp range is integrated over, 
a = 3.02 for C and 3.39 for Cu with both uncorrelated fragmentation and coalescence 
while with coalescence alone, a = 1.24 for C and 1.78 for Cu. Considering only the 
range xp > 0.1, we find a = 1.56 for C and 1.63 for Cu with fragmentation and 0.43 for 
C, 0.57 for Cu without fragmentation. The calculated a's suggest considerably harder 
Xp distributions in the more limited xp region, particularly when coalescence alone is 
considered. However, none of the results are in good agreement with obtained 
by WA89 for xp > 0.4. This is not surprising, especially since the solid curve is seen 
to be harder than the data for xp > 0.1. In any case, the coalescence contributions 
in particular, now considerably more important than for S~ production by non-strange 
hadrons, are difficult to fit to a power law since they approach zero at both xp — 
and Xp = 1 with a peak at intermediate xp, see the curves in Ref. [28]. The various 
contributions, all with somewhat different magnitudes due to the relative probabilities, 
peak at different values of a;^, complicating the situation further. 

The calculated values of a give a 0.8 for the integrated cross sections but a 0.7 
for Xp > 0.4, with and without fragmentation in the intrinsic model, rather consistent 
with the WA89 result. However, as a function of xp, coalescence alone is more consistent 
with the measurements since a~lata;i?~0, decreasing to 0.71 as a;i? 1, as expected 
from Eq. (22). 

We now check if our results improve when we inchidc flavor excitation, Eq. (14), 
shown in Fig. 4(c) and (d). Now the baseline leading twist distribution, F + X„„ , is 
harder than with fusion alone, as shown in Fig. 1(c). However, although the distribution 
is broader, it still drops six orders of magnitude over the entire xp range with pr^-^ = 1 
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GeV. Thus including flavor excitation cannot describe the data without the intrinsic 

coalescence component, as in Fig. 4(a) and (b). The total cross section is in reasonable 
agreement with that measured by WA89. Decreasing PTmin further can harden the distri- 
bution but still underestimates the data. A lower pr ■ enhances the total cross sections 

^ ^ min 

considerably so that, if Pt^^^ ~ 0.25 GeV, the cross section is overestimated by several 
orders of magnitude. The intrinsic contribution is then neghgible so that decreasing 
PTmin actually degrades the agreement with the data. Thus there is no clear evidence for 
flavor excitation. 

The trends in the A dependence are similar when excitation is included although 
the values of a obtained are suggestive of harder xp distributions than with leading- 
twist fusion alone. In particular, the excitation contribution is harder at low xp, see 
Fig. 1, causing the change in slope due to the hardening of the intrinsic distribution 
when coalescence alone is included to be less abrupt. Nonetheless, the agreement with 
the measured value of a is not significantly improved. The calculated values of ai^xp) 
are similar to those with leading-twist fusion alone but the integrated values of a are 
somewhat larger, 0.87, due to the larger leading-twist baseline. The A dependence 
also does not support flavor excitation as a significant contribution to strange hadron 
production. 

To summarize, the A dependence of S~ production by E~ is stronger because the 
intrinsic contribution with coalescence dominates the Xp distribution already dX Xp ^ 
0.1. Therefore the integrated A dependence is nearly a factor of A^^^ down relative to 
the pion and neutron data, as shown in Eqs. (22) and (23). Thus the trends of the model 
are in agreement with the WA89 data. 

Finally, we compare our intrinsic model calculations with the invariant S~ and O 
cross sections measured in S~Be interactions at 116 GeV [14]. Since intrinsic production 
is expected to be a primarily low px effect [25], we only compare to the low bin, 
< < 0.4 GeV^ [42]. The data and our calculations are shown in Fig. 5. We have 
multiplied our xp distributions by 2mTC0shy/y^ to obtain the invariant cross section. 
The invariant xp distributions are harder as a function oi xp. 

Because the initial and final states are identical for S~ production, the intrinsic 
contribution increases with xp. However, even with coalescence alone, the increase does 
not continue beyond ~ 0.4. A similar but weaker effect is seen for the Q where there 
are two strange quarks in common with the S~. Therefore we have tried to identify a 
mechanism that would increase the cross section beyond xp 0.4. One possibility is 
a "Pomeron-like" parton in the Fock state. Since the Pomeron has quantum numbers 
similar to two gluons, it can be exchanged between two projectile valence quarks. A 
\dssV) state, where V signifies the "Pomeron", would result in S~ states at high xp 
while avoiding the 5{1 — xp) delta function for the 3-particle Fock state. A S~ from 
such a configuration would have a distribution peaking at xp ^ 1. It is also possible 
to imagine a IdssssV) state from which both the S~ and Q could be produced. In this 
case, the distribution would peak away from xp ^ 1. We included "Pomeron" production 
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Figure 5: The model is compared to the 116 GeV S^Be data. In (a) and (b), the 
leading twist contribution is F while in (c) and (d), flavor excitation is also included, 
F + X„„ . The dotted curves are for leading-twist fusion alone, the dashed curves 
include uncorrelated fragmentation and coalescence, the solid curves include coalescence 
alone and the dot-dashed curves include a diffractive "Pomeron" contribution. The data 
sets have been normalized to the cross section per nucleon. The curves are normalized 
to the data at xp = 0.15. 
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from both states assuming that — P^^ ~ 5% and that i^^^ = Pj^q ~ 3.5%, giving 
these configurations large probabihty. The resuhs, shown in the dot-dashed curves in 
Fig. 5, agree relatively well with the data, especially for the 5~. The Vt data are still 
underestimated but the trend is now in the right direction. 

Of course, this is a rather artificial solution, especially when the initial and final 

states are not identical. If it is correct, it should also be included in E~p — * 
as shown in Fig. 4. However, we have checked this and found that the resulting xp 
distribution is far too hard. Therefore, the practicality of the mechanism is questionable 
and the "Pomeron" results should not be taken too seriously. 

Our calculations with flavor excitation are compared to the data in Fig. 5(c) and (d). 
The results do not improve, even when the "Pomeron" is included. Indeed, the results 
with excitation underestimate the data at high xp more than with fusion alone. There- 
fore we conclude that flavor excitation is not an effective mechanism for strange hadron 
production at low pr, in keeping with the interpretation of the excitation diagrams as 
NLO contributions to the production cross section, as discussed previously. 

We have also calculated the exponent a', see Eq. (3), for these distributions with 
Xp > 0.5. The results are negative for all the cases shown. We flnd a' ^ —0.45 
with and without uncorrelated fragmentation and —0.5 with the "Pomeron". These 
values are rather consistent with those obtained for the low selection of the S~ data. 
The corresponding values for Q production are somewhat lower, —0.41, without the 
"Pomeron" but somewhat higher, pa —0.54, with it. 



6 Conclusions 



We have compared our intrinsic calculations to S~ production by vr", n and S~ pro- 
jectiles and to S~ and fl production by S~ projectiles. We find good agreement with 
the WA89 data for leading-twist fusion and coalescence. Flavor excitation seems ex- 
cluded as a significant mechanism of low pt strange hadron production. The apparent 
difficulties with uncorrelated fragmentation seen in Ref. [28] are confirmed here. Thus 
coalescence production is the most effective mechanism for strange hadron production 
in the intrinsic model. The leading charm analysis should perhaps be revisited in light 
of this conclusion. 

The conclusions that can be reached from the S~-induced interactions at 116 GeV 
are less clear. It is possible that a "Pomeron-like" state could exist in the hadron 
wavefunction but its applicability to Q production is somewhat doubtful. Therefore 
the interpretation of these data within the intrinsic model is rather inconclusive. More 
standard studies of diffractive production in S~ Be — > E^X should be performed. 
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Appendix 



Here we give the relevant probability distributions in the intrinsic model for S~ 
and Q production used in our calculations. To more clearly distinguish between the 
probability distributions including uncorrelated fragmentation and coalescence and those 
with coalescence alone, both distributions are given. 

First, we give the distributions relevant to the WA89 measurements. We reproduce 
the S~ probability distribution from a 7r~ [28] , 

dxp 10 dxp 10 dxp 10 dxp 2 ylO dxp 7 dxp J 

I 1 ^^iL^u I 1 dP^s^a I 1 dPtJ, 1 / 1 dPZ 1 dP^,g\ 
10 dxF 10 dxF 10 dxp 2\10 dxp 12 dxp J 




1 ^-Pissd ^ ^-^isS ^ 1 / 1 dP^^ 3 dP^^ 



I " "'-^ issd 1 I ( "''^ isss I " "'-^ isss I 



10 dxp 12 dxp J 2 ylO dxp 16 dxp J 



dPf ^ l^^l^^l^^l^ .2) 
dxp 7 dxp 12 dxp 12 dxp 16 dxp 

Prom a neutron projectile, 

dPj- ^ 1 dP,f ^ 1 dPl^ ^ 1 dPl^ ^ 1 / 1 dPl^ ^ 1 dPl^\ 

dxp 10 dxp 10 rfxiT 10 (iajF 2 ylO rfxi? 13 rf.XiT' ) 

1 r/p9^' 1 r/p9^ 1 dP^^ 1 / 1 dP'^^^ 9 dP^'^'\ 

I isuu _| isud _| isdd _|_ ^ I ^ issu _|_ ^ issu \ 



10 dxp 10 (ixi? 10 dxp 2 yiO dxp 22 dxp J 

2 ^10 dxj. 22 ) 2 \10 rfxj. 2S dxp J ' ^ ^ 

c^xj. 13 dxF 22 dxp 22 dxp 28 dx^ ' ^ ^ 

The 5~ distribution from a S~ projectile is, 

dPi: 1 dppf , 1 dpp/ , 1 / 1 dP,f , 2 dP,f\ , 1 dP,^„^ , 1 dPl^ 



dxi? 10 dxp 10 (ixi? 2 \^10 dxp 6 rfxi? y 10 dxp 10 

, 1 dP^^^ ^ 1 n ^ ]_dP^ ^ 1 n ^ 3 dp;i^\ 



10 o?Xi? 2 \ 10 dxp 13 7 2 \ 10 dxp 13 dxi? 
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2 ylO d,xp 16 cix/T' y 10 d,X},- 10 rfxi? 10 dxp 10 dxi? 

1 n 1 /l^ild , 3 rfP^c \ 

2 VIO dxj. 22 dx^ ; 2 \10 dxp 22 dxi. J 
1 n 14cLg|i\ 1 }]_dl^ ^d^\ 
2\10 dxp 22 dxp J 2 \10 dxp 28 dxp J 



1 n 9 dP^ 1 n rfiJJ 12 rfigs\ 

2 \ 10 dxF 28 rfxF y 2 \ 10 dxF 32 dxp J ' 
dP^: 2 dPl^' 1 dPl?: 3 dPlS 6 rfPZF 2 c^P?£, 



dxj? 6 dxp 13 (ixi;' 13 dxp 16 cixj? 22 dxp 

I 3 c?Pig„d 14 (jPjsdd , 6 (i-F^s^ 9 ^-Pjs^ 12 o^Pjg^ 
22 22 dxF 28 dx^ 28 dxp 32 dx^ ' 



We now present the relevant probability distributions for S and Q production from 
a S~ projectile. First we give the 5~ distributions, 

dPj: ^ 1 / 1 dP^ ^ Id/^f \ ^ 1 / 1 dPPj ^ 2dP^^\ 
dxp 2 ylO d,X},^ 6 dxp J 2 \10 dxp 6 dxp J 

1 ( 1 dPif' SrfPif \ 1 / 1 dP^Z 1 (^^iuu\ 



+ 



10 dxp 7 dxp / 2 yiO dxp 13 dxp J 
1 dP"^^ 9 dP'^'^\ 1 / 1 dP'^^ 3 r/P'^'^' 



+ 



10 dxp 13 cixiT' y 2 yiO dxp 13 cixi? 
■ 1 dPll , 3 rfPiin , 1 / 1 dPl^^ , 6 rfPiir 



10 dxi? 16 dxp J 2 ylO dxp 16 dxp J 
' 1 r/P'^^ 6 dP-'^^\ 1 / 1 dP-^^ 1 r/P^"^ ^ 



10 dxi? 17 dxi;' / 2 yiO (ixi? 22 dxp / 



+ 



10 dxF 22 / 2 ^10 dxp 22 dxp ) 

"Mddd , ^ "Mddd I , / "Msuu , "Msuu \ 



+ 



10 22 dxF ) 2 ylO 28 dxp 

1 rfi^iL^d + + 1 f 1^ + 9 dP^, 



10 28 j 2\\^ dxp 28 dxjr j 

1^ ^ ^ 1 (J_d^ ^ 12 

10 dxp 31 dxi^ y 2 ylO dxp 31 dxi;' y 

1 dP^I 12 dP^S\ 



ISSS 



2 ylO (ixi? 37 dxi;' y 

^ ^ 1 rfPjf ^ 2 rfP,f ^ 3 rfP,f ^ 1 dP,I^ ^ 2 ^ 3 rfP,^,g 

dxp 6 (ixi? 6 rfxi? 7 rfxi? 13 rfxj? 13 dxp 13 rfxi? 

, "J isu , " "Msd , " iss _|_ iuuu _|_ ^ iuud _|_ "Mudd 



16 (ixi? 16 dxp 17 (ixi? 22 dxp 22 dxp 22 dxp 
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22 28 rfx^ 28 dxp 28 cix^ 

, 6 (iPis^j 12rfi^g^ 12 (iPjg2 



31 dxp 31 (ixf 37 
Finally, we give the Q distribution from a S~ projectile, 

dPE' ^ 1 rfPjf ^ 1 dPt/ ^ 1 / 1 rfPjf ^ IciPif \ ^ 1 dP^,^ ^ 1 ciP,^,^ 
dxp 10 (ixi? 10 rfa^i? 2 \10 dxp 7 dxp J 10 dxp 10 (ia;^ 

1 ^-^ild _|_ 1 dPisd\ 





16 dxp J 2 ylO dxp 16 dxp J 



10 dxp 17 / 10 dxp 10 10 dxiT 10 dxp 

10 dxp 28 dx^^ ; 2 \10 dxF 28 dxj. j 
1 dP^^^^ ^ 1 dP^g^^ \ ^ 1^ /j^rfij^ ^ 3 (iPjg£^ 



2 yiO cixiT' 28 dxp J 2 ylO dxp 31 dxp J 

2\10 dxp 31 dxp ) 2 VIO dxp 37 / ' 

dxp 7 dxp 16 (ixiT' 16 dxp 17 lixi? 28 dxp 28 

, 1 '^-^isdd _|_ 3 c^Pjs^ _|_ 3 (iPjg^ ^ 10 (i-F^g^ ^^^^ 



28 cixi;' 31 (ixi? 31 dxp 37 (ia;^ 
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